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A thin outer ice shell on Jupiter’s large moon Europa would imply 
easy exchange between the surface and any organic or biotic 
material in its putative subsurface ocean 1-4 . The thickness of the 
outer ice shell is poorly constrained, however, with model- 
dependent estimates ranging from a few kilometres 5 6 to ten or 
more kilometres 7 . Here I present measurements of depths of 
impact craters on Europa, Ganymede and Callisto that reveal two 
anomalous transitions in crater shape with diameter. The first 
transition is probably related to temperature-dependent ductility 
of the crust at shallow depths (7-8 km on Europa). The second 
transition is attributed to the influence of subsurface oceans on 
all three satellites 3,8,9 , which constrains Europa’s icy shell to be at 
least 19 km thick. The icy lithospheres of Ganymede and Callisto 
are equally ice-rich, but Europa’s icy shell has a thermal structure 
about 0.25-D.5 times the thicknesses of Ganymede’s or Callisto’s 
l shells, depending on epoch. The appearances of the craters on 
Europa are inconsistent with thin-ice-shell models 1 and indicate 
that exchange of oceanic and surface material could be difficult. 

The thickness of Europa’s ice shell controls whether exchange 
between Europa’s surface and its putative water ocean (including 
exposure of possible organic materials to sunlight 1 ) is possible, and 
in what manner. A thick shell (>10 km thick) would be susceptible 
to solid-state convection and overturn 10 , whereas a thin ice shell 
(only a few kilometres thick) would be vulnerable to crack-through 
and melt-through from below 1,5 . Efforts to constrain the thickness 
of Europa’s icy shell through geologic 3-7 and geophysical 3,4 investi- 
gations have thus far failed to produce a consensus because of 
controversy over the origin of geologic structures 1,2 and weak model 
constraints. Impact craters (Fig. 1) sample planetary lithospheres 
over a range of depths and offer direct probes into Europa’s icy shell. 
Crater morphology is influenced by planetary surface gravity" and 
lithospheric properties 12,13 . Surface gravity on the three icy galilean 
satellites and the Moon are similar (Europa, 132cms - "; Ganymede, 
143 cm s” 2 ; Callisto, 124 cm s -2 ; the Moon, 162 cm s -2 ), so differ- 
ences in crater morphology can be attributed to compositional or 
rheological differences. 

Moore et u/. 14,1 \ in their initial survey of Galileo observations of 
crater morphologies on Europa, found that the largest craters (Fig. 
1) are morphologically distinct from and shallower than similar- 
sized craters on Ganymede and Callisto (Fig. 1). They concluded 
that these craters may have excavated into either a global subsurface 
ocean or a ductile zone in the lower shell, but did not conduct a 
i systematic analysis of crater dimensions. Recent attempts to model 
I impact melt production suggest that the ice shell may be at least 3- 
| 4 km thick 16 , but this does not constitute a definitive test of the ‘thin 
crust, thick crust’ debate. 

Depth/diameter (d/D) data for fresh, undegraded, craters on the 
1 Moon and other planets follow well-defined trends and feature two 
breaks in slope 11,17 . These breaks, or transitions, correlate with 
morphologic transitions from simple bowl-shaped to complex 
morphologies (featuring slump terraces and uplifted central 
peaks) and from complex craters to ringed basins (see, for example, 
ref. 18). For comparison, I have catalogued all craters on Ganymede 
and Callisto that are over 30 km across, and over 1 km across on 
Europa. I have also measured the depths of fresh (unmodified and 


unrelaxed) craters from 100 m to 170 km in diameter on the icy 
galilean satellites (Fig. 2) using Voyager and Galileo images. Three 
transitions in crater shape have been identified on each satellite, but 
two of these are anomalous and occur at different diameters on 
Europa than on Ganymede and Callisto. 

Transition I (from simple to complex morphology) is similar on 
all three satellites (Fig. 2) and can be associated with regular crater 
modification processes 18 . Transition II on these satellites correlates 
with anomalous changes in complex crater dimensions and mor- 
phology. On Ganymede and Callisto, this occurs at D ~ 26 km, 
where crater depths begin to decrease slightly (d ~ 1.1 ± 0.3 km) 
and are shallower than predicted from extrapolation of the d/D 
curve for complex craters (Fig. 2a, b). This change coincides with 
the morphologic transition from central peak craters to central pit 
and dome craters 19 (Fig. 1). Transition II on Europa occurs at 
D - 8 km, where complex craters on Europa diverge from their 
galilean siblings and become increasingly shallower than similar- 
sized complex craters on Ganymede and Callisto (Fig. 2c). Central 
peaks still occur but crater rim and peak morphology becomes 
progressively disrupted with diameter (Fig. 1). 

Transition III is associated with a sharp reduction in crater depths 
and the development of anomalous impact morphologies. On 
Ganymede and Callisto, many large central dome craters have 
poorly developed or non-existent rims (Fig. 1). Floors of these 
‘anomalous dome craters’ (Fig. 1) are similar to or elevated a few 
100 m above surrounding plains. Rim-to-floor depths are negligible 
to about 500 m (Fig. 2b). Anomalous dome morphologies in older 
craters (based on superposed crater densities) occur at D > 60 km, 
indicating a Transition III diameter as low as 60 km in post-bright- 
terrain formation times. For younger central dome craters only 
those larger than about 150 km in diameter have anomalous 
morphologies, suggesting that Transition III diameters on Gany- 
mede and Callisto have increased to about 150 km at present. 

Transition III on Europa occurs at D = 30km (Fig. 2c) and 
correlates with an abrupt change from modified central peak to 
multiring morphologies between 27 and 33 km crater diameter. The 
two largest impact structures, Callanish and Tyre (D ** 33 and 
41 km, respectively, as determined by scaling from mappable ejecta 
deposits 20 ) are characterized by several ( >5) concentric ridges and 



Figure 1 Impact crater landforms on the icy galilean satellites, representing the major 
mo^phologic types. In order of increasing diameter from left to right crater types are, for 
Ganymede and Callisto (top row): central peak craters (D = 18 km), central pit craters 
(D ---- 30 km), central dome craters (Enkidu, D = 121 km) and anomalous dome craters 
(D == 138 km). For Europa (bottom row): central peak craters (0 = 8 km), modified 
central pit craters (Grainne, D — 14 km) r anomalous central peak craters (Pwyll, 

D ■ 27 km) r and multiring basins (Tyre, 0=41 km). Images have been scaled so that 
craters appear to have the same diameter. Scale bars are 30 km for Ganymede and 
Callisto (top row), and 10 km for Europa (bottom row). Images for Pwyll and Tyre and the 
unnamed Ganymede basin have been colour-coded to represent topography (red, high; 
blue, low; colour scale bars represent relative elevation: 0 to 500 m for Europa craters and 
0 to 2,000 m for Ganymede). Note the progressively more complex central peak 
mc r phology with increasing crater diameter, The nominal rim location for Tyre lies within 
the nnermost of the prominent topographic ridges. The Sun is from the left in all images. 
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graben 13 , but no central structures, depressions or rim scarps can be 
identified within the nominal crater (Fig. 1). Maximum effective 
crater depths, inferred from isolated ridges near the rim location 
and subtle topography across the floor, are between 50 and 100 m 
(Figs 1, 2c). These structures differ sharply from terrestrial planet 
multi ring basins which have only 3-5 ring scarps and topography 
greater than 4 km (for example, refs 13, 17). 

The rollover and precipitous drop in crater depths with diameter 
for the galilean satellites (Fig. 2) is unprecedented among the 
planets. That this is not due to long-term post-impact flattening 


! 




Figure 2 Depth/diameter measurements for fresh impact craters on the icy galilean 
satellites. Data are plotted separately for Callisto (a), Ganymede (b) and Europa (c). The 
heavy dark line is for lunar craters 11 ,7 . The thin lines are least-squares fits through data. 
Craters plotted are simple craters (solid dots), complex (central peak} craters (open 
circles), central pit and central dome craters (crosses), and anomalous dome and 
multiring basins (error bars). For Europa, the modified central peak craters Mannann'an 
and Pwyll are plotted as circles with crosses. The dashed line in c is the least-squares fit 
through complex craters on Ganymede for comparison with complex craters on Europa. 
Note the sparser data set for Callisto. The error bar for the largest central dome crater on 
Callisto is typical for all plotted data. 


by creep (that is, viscous relaxation) is suggested by the youth of the 
craters studied (some have well-preserved bright rays) and by their 
misshapen or absent central structures and rims. Relaxation and 
volcanism may flatten or bury topography but are unlikely structu- 
rally to disrupt rims and central peaks (see ref. 21). Rugged relief on 
crater floors (Fig. 1) is also inconsistent with the idea that craters 
formed within a giant pool of impact melt 16 . Thus, differences in 
internal structure between the icy satellites and the Moon, and 
between the icy galilean satellites, are the most plausible explanation 
for the unusual shapes of larger craters on these satellites. Warm ice 
probably at shallow depths within these satellites is notoriously 
weaker than common silicates and the very cold brittle ice near their 
surfaces (see ref. 22). 1 propose that Transitions II and III reveal at 
least two such temperature-induced rheologic transitions or phase 
changes with depth on all three icy galilean satellites, but that these 
transitions are much shallower on Europa. Inherently weaker ice at 
greater depths would be progressively less able to support the 
topography associated with larger transient craters several kilo- 
metres deep 23 , leading to enhanced structural collapse. 

Impact into Theologically layered targets is imprecisely under- 
stood. Studies of terrestrial and lunar craters 18,24 and recent model- 
ling 23 of crater floor rebound suggest that modification from simple 
to complex morphology (by floor rebound and rim collapse) 
involves a quasi-hemispherical zone around the transient (pre- 
modification) crater 18 and that the modification process and final 
crater shape begins to be unduly influenced when an anomalously 
weak layer is roughly 1 to 1.3 times as deep as the transient crater 
width. The McKinnon and Schenk 23 scaling relationship for icy 
satellite craters is used to estimate transient crater dimensions (to 
within approximately 10%) from the observed crater diameters. On 
Ganymede and Callisto, Transition II craters ( D ~ 26km) scale to 
transient craters about 16.5 km wide for a depth to the first 
anomalous layer of —16-22 km. On Europa, Transition II craters 
{D ~ 8 km) scale to transient craters about 6km wide, for a layer 
depth of 7-8 km. On Ganymede and Callisto, Transition III craters 
(D ~ 150 km) scale to transient craters about 80km wide, for a 
depth to the second anomalous layer of 80-105 km (the older 
Transition III diameter of — 60 km translates to transient diameters 
of around 35 km, and a layer depth of 35-45 km). On Europa, 
Transition III craters (D ~ 30km) scale to transient craters about 
19 km wide, and a depth to the second transition of 19-25 km. 

Transition II is interpreted to represent a temperature-dependent 
rheologic change with depth. The preservation of prominent (albeit 
modified) central peaks in complex craters not yet at Transition II is 
consistent with solid-state material at the base of the transient 
crater 14,16 . Rounded central domes, modified central peaks and 
degraded rim topography in craters larger than those at Transition 
II (Fig. 1 ) are all consistent with coherent but very warm ductile ice 
in the lower icc shell (see ref. 19). 

Galileo tracking data indicate that the interior of Ganymede is 
probably strongly differentiated with an ice-rich lithosphere and 
mantle, whereas Callisto is not fully differentiated and may or may 
not have an ice-rich outer lithosphere 26,27 . The indistinguishable 
dID statistics for craters on Ganymede and Callisto (Fig. 2) indicate 
that internal composition and thermal structure of these two 
satellites are similar, and that Callisto also has an ice- rich outer 
shell at least several tens of kilometres thick. 

Galileo magnetometer observations are consistent with internal 
oceans currently within both Ganymede and Callisto at depths of 
150-200 km hut possibly shallower 8,9 . If Transition III is a detection 
of oceans on Ganymede and Callisto, they must be at least 80 km 
deep at present. Changes in morphology and transitions on Gany- 
mede and Callisto as described here may record changes in litho- 
spheric thermal properties 28 and possibly ocean depths over time. 

Transition III to multiring morphologies on Europa may indicate 
the beginning of the influence of a global ocean 13 . The unusual 
multiring morphologies on icy satellites are probably related to 
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impact into a thin brittle lithosphere overlying a fluid or a ductile 
layer that behaves as a fluid during crater collapse 10 . The apparent 
abruptness of Transition III on Europa ( Figs 1,2) is consistent with a 
sharp phase change at depth, such as an ice-liquid interface, but 
may also represent a rapid transition to the hot base of the ice shell 
where the ice is near melting and very ductile. It seems that the base 
of the icy shell is not much deeper, but the inferred Transition III 
depth at least indicates a minimum ice shell thickness of around 1 9- 
25 km. Europa’s known impact craters are unlikely to excavate 
ocean material directly to the surface, however (unless the material 
is trapped in isolated shallow reservoirs). The largest crater, Tyre, 
excavates down to only about 3 km (~30-40% of the transient 
crater depth 1 *), or around 15% of the minimum shell thickness. 

Systematic mapping of impact crater shape and morphology 
provides a robust (and free) means of sampling icy satellite 
interiors. Ice rheology is dependent on temperature 22 and these 
transition diameters provide new constraints on the comparative 
thermal structure of these satellites. Transitions II and III are 2—4 
times shallower on Europa than on Ganymede or Callisto, indicat- 
ing that the rheologic structure of Europa s outer shell is similarly 
thinner than on Ganymede and Callisto. As heat flow scales roughly 
linearly with the thickness of the stagnant non-convecting lid 10 
(possibly represented by Transition II), heat flow should to first 
1 order also be 2-4 times higher on Europa, depending on epoch. 

: Transition III constrains Europas ice shell to be at least 19-25 km 
thick, consistent with that required for convection to proceed 
within the ice shell 10 and favouring diapirism of the lower shell 
for the origin of ovoid features and chaos terrains (refs 7, 29 and P.S. 
and R. Pappalardo, manuscript in preparation). The minimum shell 
thickness strongly supports thick ice shell models in general, and 
associated interpretations of geologic features 2 ' 29 ; it may also 
indicate the depth to Europas ocean, or at least the beginning of 
the hot basal layer of the floating ice shell. □ 

Methods 

Crater depths have been measured from Galileo and Voyager data using three techniques: 
stereo digital elevation models (OEMs), photoclinometry OEMs, and shadow length 
measurements. Stereo DEMs resolve only the large craters Pwyll, Mannann’an, and t.ilix 
(diameters 19-27 km) on Europa, and several anomalous dome craters on Ganymede. 
Photoclinometry in two horizontal dimensions can be used with single low-sun images to 
map topography from relative brightness. My photoclinometry technique includes the use 
of low-phase angle images to model local albedo, thus reducing or eliminating one of the 
major systematic sources of error in photoclinometry. This technique is used here 
primarily to confirm the stereo and shadow measurements, but is the primary source of 
topographic data for the multiring basins Callanish and Tyre. Photoclinometry was also 
used to map topography across additional anomalous dome craters on Ganymede, 
j supporting the conclusion that they have raised floors based on limited stereo data 
j described above. 1 also use low- resolution controlled stereo DEMs to control the long- 
wavelength component of high- resolution photoclinometry DEMs. Phis was especially 
; useful for Pwyll (Fig. 1 ). Height measurements based on triangulation of shadow lengths 
are used for all craters on Europa, Ganymede and Callisto smaller than ~ 15 km across. 
Depth measurement errors include systematic errors due to technique and those due to 
variations in topography along the rim crest. Systematic errors rarely exceed 10%. 
Variations in rim height approach about 100 m for large craters such as Pwyll. These data 
supercede the depth/diameter statistics for Ganymede and (Callisto based on lower- 
resolution Voyager data' 2,1 *; no Voyager-based measurements were possible for craters on 
Europa. 
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